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UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



Aux* 
LIBRARY ROUTINE E6-221 



Integration over a Single Interval (Tabulated Values ) 

Closed subroutine 

32 (see Note) 

0, 1 

(Quadrature (see Method) 

3.7 + l-6s msec. 

This routine computes an approximation to the integral 



f (x) dx/w 



when f (x) is tabulated for the (2s +2) values, x 




sw, a - (s-l)w, - - -, a^ - w. a 



a + w. - - «-, 



a + sw, a + (s + l)w. 
P P 



P * P' P 
The order of the integration 



process and therefore its accuracy, for given interval w 
and function f (x), is controlled by the integration 
parameter s . The location of f , = f (a + kw ) in the memory 
must be (p + k) where p is arbitrary and | f | < l/2 . 
The table of f (x) must extend at least s intervals to the 
left and to the right of any interval over which f (x) is 
to be integrated. 

The two parameters, p and s, appear in the entry to the 
subroutine as follows: 





50 p F 


q. 


50 q F 




26 — 


q + 1 


00 s F 



Control is returned to the left hand side of (q + 2) with 



the result. / 

Va 



P 



+ w 



f (x) dx/w, in A and also in location 



0. The parameters a and w need not be specified since the 
result is unchanged by a linear transformation of the variable 
x. 



_2- 



Since the subroutine permits efficient numerical integration 
with evenly spaced ordinates, it is specially useful when 
the interval w is large because of prohibitive cost of 
computing ordinates . Also* because s appears in the entry, 
the subroutine can be used for calculations with several 
values of s chosen by another routine, as required for 
example when integrations must be repeated until the error 
is reduced below a pre-set value. 

Such an application might occur, for example , during the 
calculation of an indefinite integral, 

NOTE The 10 constants needed for calculations with s = 0, 1, 2, 

3 are stored at 22L, 23L «— 3-U* in the subroutine. If a 
value of s = s 1 > 3 is to be specified in the entry then 
the corresponding (s 1 +1) constants must be placed at relative 
positions (u n + 22) L, (u., + 23)L,- - - (u., + 22 + s )L where 
vl = s,(s 1 + l)/2o If the largest value of s needed , 
say s = s < 3, then the subroutine can be shortened by 
omitting constants at the end of the subroutine, 

METHOD The numerical integration formula* used is 



/ 



a + w s 



, & f (x) dx/w = E A ±s (f^ + f ± + ± ) + e(s,v) = g g + e(s,w) 

The quadrature error term e(s,w) corresponds to approximating the 
function f (x) by a polynomial of degree (2s + l) and can be 
estimated from 

|c(b,v)| < ( W /2) 2S +2 | f (2s + 2) (a + | )| , |J | <» 

The coefficients A. can be generated from the equation 
is 

g s -8 S -l= 2P s* B2Sf l/2 

2s 
where u5 f , is the mean centered difference of degree 

2s at x = a + w/2 and 



■^Integration of the Newton^Bessel interpolation formula (Whittaker and Robinson 
(2nd Ed.), p, ko) gives a quadrature formula that can be expressed in terms of mean 
central differences (loc. cit., p. lVf). In the formula used by the subroutine, the 
central differences have been replaced by appropriately weighted sums of ordinates . 



EXAMPLE 



-3- 



P = 
s 



[2 2S + 1 (2s)i] - 1 ACz 2 - l)(z 2 - 9)— [z 2 - (2s - I) 2 ] dz 
^0 



Note that A = P . 

ss s - 

The error in calculating / sin x dx/w with this 

subroutine, for a = 55 Jt/l80, v = 15 ot/l80 and s = 3> is 
6.6x10 « The error estimated "by the above formula is 
|e(3, 15 it /l80)| = 7.6x10 . 
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4l IF 
K5 F 








1 


42 4L 
46 11L 






— Set addresses 


2 


LO 20L 
10 20F 








3 


42 11L 

F5 4l 








4 


42 l8L 




Plant link 




L5 F 








5 


42 IF 
4l F 








6 


50 IF 
75 IF 






— Set addresses 


7 


S5 F 
L4 IF 








8 


10 IF 
L4 21L 








9 


42 12L 
L4 IF 








10 


42 19L 




Set end constant 




36 18L 




End test 


11 


L5 F 
L4 f 








12 


40 IF 
50 L 






r f .i + f i + i to1 


13 


TJ IF 
l4 f 






" 2 A ls (f_. + t . +z ) to 


14 


40 F 
L5 11L 






r 


15 


L4 20L 
40 ITT, 






— Change addresses 


16 


F5 12L 
42 12L 








IT 


FO 19L 
22 10L 






— Prepare for end test 



e6 



LOCATION 


OEDER 
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PAGE 2 


18 


L5 F 
26 F 




Sum to A 
Link 




19 


to IF 
50 L 










20 


LL 4095F 
00 IF 






— Constants 




21 


80 F 
00 22L 










22 


to F 










00 0000 0000 OOOOJ 


^00 




23 


to F 










oo oki6 6666 666yj 


A 01 




2k 


NO F 










00 4583 3333 3333J 


A ll 




25 


to F 










00 0569 khkk hkkhj 


^02 




26 


NO F 










00 hj^h 1666 6667J 


A^ 




27 


00 F 








00 0076 3888 8889J 


a 
22 




28 


to F 










00 0648 3961 6to2J 


^3 




29 


NO F 










00 4212 0535 71^3J 


*13 




30 


00 F 










00 0155 I 


to6 0847J 


A 23 




31 


NO F 










00 4984 2 


096 5608J 


A 

33 





e6 



